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J.  Awerbuch  and  S.R.  Bodner 

Abstract 


A  mathematical  model  has  been  developed  to  describe  the  mechanism  of 
normal  perforation  of  projectiles  in  metallic  targets.  The  perforation  process 
is  considered  to  be  divided  into  three  interconnected  stages.  The  analysis 
accounts  for  an  effective  mass  of  the  bullet  due  to  part  of  the  target  material 
moving  with  the  bullet,  the  deformation  of  the  bullet  during  penetration,  and 
the  increased  strength  of  the  target  material  at  high  rates  of  loading.  The 
analysis  enables  the  residual  velocity  to  be  calculated  as  a  function  of  the 
target  thickness  and  its  mechanical  and  physical  properties,  and  of  the  mass, 
geometry  and  impact  velocity  of  the  projectile.  The  geometry  of  the  cavity, 
i.e.  entrance  and  exit  diameters  and  plug  thickness,  are  factors  in  the 
analysis  and  are  empirical  quantities.  The  present  theory  can  also  predict 
the  force-time  curve  ana  the  contact  time  for  the  perforation  process. 
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Notation 


projection  of  nose  of  projectile  on  the  target  plate 
cross  sectional  area  of  the  cavity  in  the  first  stage 


-  cross  sectional  area  of  the  cavity  in  the  second  stage 
A  -  cylindrical  surface  area  of  plug 
b  -  plug  length 
D 
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D. 


diameter  of  the  cavity  in  the  first  stage  of  the  perforation  process 
(considered  equal  to  entrance  diameter) 

diameter  of  the  cavity  in  the  second  stage  of  the  perforation  process 
(diameter  of  plug) 


C3  -  exit  diameter 

e  -  radial  width  of  shear  zone  of  the  target  plate 
F  -  resultant  force  on  effective  mass  of  projectile 
•  inertial  force 
Fc  -  compressive  force 
Fs  -  shearing  force 

Fj,  ^2 <  "  total  forces  acting  on  the  combined  projectile  and  effective 

1  added  mass  during  the  different  stages  of  the  perforation  process 

Fx2‘  ^c2*  1 s2  *  inertial>  compressive  and  shearing  forces  acting  on  the 
projectile  during  the  second  stage 

li  -  thickness  of  target  plate 

K  -  numerical  constant  depending  on  the  shape  of  the  projectile  nose 
m  -  instantaneous  mass  of  projectile 
ms  -  original  mass  of  projectile 
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Notation  (cont'd) 

-  projectile's  mass  at  the  end  of  the  first  stage 

m2  -  projectile's  mass  at  the  end  of  the  second  stage- 

t  -  time 

t, ,  t2,  tj  -  duration  times  for  the  different  stages  of  perforation 
t^.  -  total  time  for  the  perforation  process 
V  -  instantaneous  velocity 

-  impact  velocity 
Vj.  -  final  velocity 

Vj,  V2,  Vj  -  projectile's  velocity  during  the  different  stages  of  the 
perforation  process 

*  velocity  at  end  of  second  stage 

9^  -  velocity  at  end  of  third  stage,  equal  to  V^. 

penetration  depth  of  the  projectile  and  effective  added  mass 

semi  apex  of  conical  nose  of  projectile 

dynamic  ultimate  shear  strain 

shear  strain  rate 

p  -  density  of  target  material 

U  -  coefficient  of  viscosity  for  shearing  deformation 
c>c  -  dynamic  ultimate  compressive  stress 
t  -  dynamic  ultimate  shear  stress 

Q  -  displacement  of  combined  projectile  and  plug  during  third  stage 
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Introduction 


Perforation  of  a  target  plate  due  to  the  impact  of  a  projectile  may 
occur  by  a  number  of  mechanisms  such  as  petal  formation  (or  dishing),  ductile 
hole  enlargement,  plug  formation,  and  the  fragmentation  (scabbing)  of  the 
target  material  (as  shown  schematically  by  Goldsmith  [1]).  Various  theories 
have  been  proposed  to  explain  the  resistance  of  metallic  plates  to  projectile 
penetration.  Due  to  the  complexity  of  the  problem,  the  suggested  analytical 
models  are  generally  simplified  by  some  basic  assumptions  anu  appioximations . 

The  two  main  approaches  that  have  been  used  to  analyze  this  problem  are 
those  of  energy  balance  and  of  conservation  of  momentum.  The  energy  balance 
method  was  applied  by  Taylor  [2]  who  studied  the  enlargement  of  a  circular 
hole  by  a  conical  head  projectile  perforating  a  thin  plate  and  derived  an 
expression  for  the  total  work  required  for  plastic  deformation.  Thomson  [3], 
also  using  the  energy  method,  derived  equations  for  the  energy  dissipation 
due  to  plastic  deformation,  heating,  and  inertial  resistance  of  the  target 
material.  A  similar  approach  was  proposed  by  Brown  [4]  to  evaluate  the 
energy  dissipated  during  the  process  of  bullet  containment  in  thin  plates, 

A  different  approach  for  the  case  of  perforation  of  thin  plates  was 
proposed  by  Zaid  6  Paul  [5,6].  Their  method  is  based  on  momentum  balance 
for  the  target-projectile  system  which  requires  that  the  terminal  shape  of 
the  perforated  plate  be  specified.  A  similar  procedure  was  used  by 
Nishiwaki  [7]  who  proposed  a  theory  for  the  perforation  of  thick  plates 
based  upon  data  derived  from  static  tests. 

Investigations  on  plug  formation  during  perforation  were  made,  among 
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others,  by  Recht  6  Ipson  (81  who  dealt  with  the  case  of  high  velocity  im¬ 
pact.  Recht  5  Ipson  (8)  developed  an  eneigy  analysis  for  the  case  of  the 
plug  mode  of  failure  from  which  the  residual  velocity  could  he  calculated 
provided  the  minimum  perforation  velocity  is  known. 

Most  proposed  analyses  are  restricted  to  tho  case  of  high  velocity 
impact  in  order  to  justify  a  number  of  assumptions  such  as  a  constant  velocity 
during  the  perforation  of  thin  plates,  the  absence  of  plastic  deformation 
beyond  the  immediate  rone  surrounding  the  hole,  and  a  constant  pressure  on 
the  projectile.  These  analyses  are  also  restricted  to  the  case  for  which 
the  projectile  is  not  deformed  during  the  perforation  process  and  are  gene¬ 
rally  based  upon  only  one  of  the  possible  mechanisms  of  perforation.  Actual 
perforation  of  a  target  plate,  however,  may  occur  by  a  combination  of  two  or 
more  mechanisms.  For  example,  the  thickness  of  the  sheared  plug  is  generally 
smaller  than  the  target  thickness  and  the  plugging  process  comir.nnces  only 
after  the  projectile  is  embedded  some  distance  in  the  target  plate.  For  this 
reason  an  upproach  based  on  a  single  deformation  mechanism  would  not  be 
applicable  for  the  case  of  projectile  impact  at  ordnance  velocities.  In 
fact,  the  residual  velocities  derived  from  those  theories  based  on  a  single 
mechanism  are  generally  higher  than  the  experimental  results,  e.g.  (9). 

The  present  investigation  is  an  analysis  of  the  perforation  process 
considering  various  deformation  mechanisms  to  be  acting  at  different  stages 
of  the  process  In  this  manner  the  various  types  of  deformations  could  be 
considered  m  on  overall  manner  and  the  analysis  would  be  more  representative 
of  the  actual  circumstances.  A  preliminary  study  of  this  nature  was  perform¬ 
ed  by  Awerbuch  (10]  and  expanded  upon  by  Goldsmith  and  Finmg.  n  (11]  .  The 
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present  papei  is  a  further  development  of  those  investigations  and  resolves 
a  number  of  limitations.  The  perforation  process  is  considered  to  consist 
of  three  interconnected  stages.  The  present  analysis  enables  fairly 
accurate  predictions  of  post  perforation  velocities,  contact  times,  and  force- 
time  histories.  The  analysis  still  relies  on  a  few  empirical  quantities 
which  can  be  determined  from  p  small  number  of  tests.  Once  these  are  deter¬ 
mined  for  a  given  projectile  and  target  material,  predictions  can  be  obtain¬ 
ed  over  a  wide  range  of  projectile  velocities  and  target  thicknesses. 
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Preliminary  Discussion 

A  relatively  simple  model  to  describe  the  mechausm  of  penetration  and 
perforation  of  projectiles  in  metallic  plates  was  presented  in  an  earlier 
paper  (10).  In  that  formulation,  t lie  perforation  process  was  divided  into 
two  stages.  The  first  was  the  compressive  stage  in  which  the  forces  acting 
on  the  projectile  were  an  ineitial  force  and  a  compressive  force,  and  the 
second  stage  was  that  of  plug  formation  and  ejection. 

The  inertial  force  in  the  first  stage  is  due  to  the  acceleration  of  the 
mass  of  the  target  material  in  contact  with  the  projectile  in  the  direction 
of  motion.  The  expression  for  this  force  component  is  obtained  by  equating 
the  worK.  done  by  the  inertial  force  acting  on  the  projectile  to  the  change 
of  kinetic  energy  of  the  displaced  target  material.  The  compressive  force 
also  acting  on  the  projectile  is  due  to  the  compressive  strength  of  the  tar¬ 
get  material  in  contact  with  -.he  projectile.  Another  basic  assumption  for 
this  stage  is  that  mass  from  the  target  material  is  added  to  the  proj^.i-ile 
during  the  penetration  process. 

The  second  and  final  stage,  according  to  the  preliminary  analysis, 
starts  when  the  ejected  plug  is  set  into  motion  as  a  rigid  body.  In  this 
stage,  the  inertial  fuice  and  tiie  compressive  force  do  not  act  on  the  pro¬ 
jectile  and  the  mass  of  the  projectile  does  not  change.  The  only  force, 
therefore,  during  this  stage  is  that  due  to  the  shearing  of  the  ejected 
plug  from  the  target  plate  which  wus  assumed  to  be  constant. 

Final  velocities  of  projectiles  computed  according  to  this  preliminary 
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model  were  compared  to  experimental  final  velocities  and  reasonably  good 
agreement  was  obtained  when  the  mechanical  strength  properties  were  raised 
in  an  arbitrary  manner  to  account  for  the  high  rate  of  straining.  The 
comparison  was  carried  out  for  the  case  of  0.22  inch  caliber  lead  bullet 
having  a  muzzle  velocity  of  400  m/sec  and  target  plates  of  commercially  pure 
aluminum,  aluminum  alloy,  and  mild  steel  of  1-6  mm  thickness. 

The  results  of  an  extensive  series  of  ballistic  experiments  have  been 
reported  by  Goldsmith  and  Finnigan  (11  j  Hard  steel  spheres  of  0.123  to 
0.S  inch  diameter  impacted  and  perforated  0.0S  to  0.25  inch  thick  2024 
aluminum  and  SAE  1020  and  4130  steel  alloy  plates  at  impact  velocities  of 
500  to  8800  ft/sec.  Comparisons  were  made  of  the  experimental  final  velocities 
to  Those  calculated  on  the  basis  of  the  preliminary  model  UO]  ana  to  a 
slightly  modified  ferm  of  it  developed  in  [llj.  The  principal  modification 
was  that  the  shearing  force  in  the  second  stage  of  the  perforation  process 
was  taken  to  be  proport  1  to  the  length  of  the  plug  still  in  contact 

with  the  target  plate  r  r  than  a  constant.  Comparison  of  the  final  velo¬ 

city  obtained  with  this  modification  to  that  obtained  experimentally  showed 
slightly  better  agreement. 

The  experimental  results  of  [llj  and  those  of  an  extensive  experimental 
program  recently  lOnoucted  (iSj  have  indicated  that  the  preliminary  model 
developed  in  (10J  is  incomplete  on  some  important  points.  These  are: 

(1)  The  final  velocities  of  projectile  computed  from  the  model  are 
not  in  good  agreement  with  experimental  final  velocities  for  the  cases  of 
high  velocity  piojcctiies  and  ru:  perforation  of  thick  plates. 
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C^)  The  preliminary  model  does  not  predict  an  important  experi¬ 
mental  observation,  namely,  that  the  difference  between  the  initial  and 
final  velocities  decreases  for  initial  velocities  slightly  greater  than 
the  ballistic  limit. 

(3)  The  force-time  curve  obtained  from  this  model  is  not  fully  real¬ 
istic. 

In  order  to  overcome  these  limitations  of  the  preliminary  model,  a  more 
detailed  analysis  has  been  developed  in  which  the  perforation  process  was 
considered  to  consist  of  three  separate  but  interconnected  stages.  The 
observations  and  results  of  an  extensive  experimental  program  served  to 
motivate  certain  assumptions  of  the  analysis.  Those  experimental  results 
and  comparisons  with  predictions  bi.sed  on  the  present  analysis  are  presented 
in  an  associated  paper  [15] . 
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Analysis 

One  principal  mechanism  of  the  plate  perforation  process  is  plug 
formation  and  ejection.  When  the  thickness  of  the  sheared  plug,  b,  is  equal 
to  ter c,  there  is  no  punching  and  the  mechanism  of  failure  is  considered  to 
be  the  ductile  type.  When  the  plug's  thickness  is  equal  to  the  thickness 
of  the  target  plate,  h,  perforation  is  completely  by  plug  formation.  However, 
the  mechanism  of  perforation  is  usually  a  combination  of  the  two  processes 
with  a  transition  stage  between  them.  The  ratio  of  the  plug  to  plate  thick¬ 
nesses  b/h  depends  primarily  on  the  mechanical  and  physical  properties  of 
the  projectile  and  target  materials. 

In  the  first  stage  of  penetration,  hig.  la,  shearing  does  not  occur  so 
this  stage  is  identical  to  that  described  in  the  preliminary  analysis  [10]. 
The  only  forces  acting  on  the  projectile  are  the  inertial  force  and  the 
compressive  force 

The  second  stage  of  penetration.  Fig.  lb,  rs  the  onset  of  shearing  of 
a  plug  from  the  target  plate.  In  this  stage  of  incipient  plugging,  three 
forces  are  considered  to  act  on  the  projectile:  an  inertial  force  and  a 
compressive  force  (as  formulated  in  the  first  stege) ,  and  a  shearing  force. 
The  shearing  force  is  due  to  the  motion  relative  to  the  target  plate  of 
target  material  which  is  accelerated  by  the  projectile  during  this  stage. 

The  change  of  the  effective  mass  of  the  projectile,  due  to  the  addition  of 
target  material  moving  with  it,  is  considered  in  this  stage  as  well  as  in 
the  first  stage  This  stage  ends  when  the  plug  is  completely  joined  to  the 
projectile  and  both  arc  moving  at  the  same  velocity.  Fig.  lc. 
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The  third  stage.  Fig.  Id,  starts  when  the  ejected  plug  and  the  pro¬ 
jectile  are  moving  together  as  a  rigid  body.  The  only  force  during  this 
stage  is  the  shearing  force  wlach  acts  on  the  plug's  circumference  and  along 
its  whole  length.  The  viscous  nature  of  this  shearing  force  at  high  rates 
of  straining  is  consideied. 

Since  the  time  of  contact  between  the  projectile  and  nie  target  is  very 
short  (the  time  duration  is  about  10-30  usee  for  the  case  of  a  0.22  inch 
caliber  lead  bullet  moving  at  a  velocity  of  400  m/sec) ,  the  heat  generated 
at  the  projectile  -  target  plate  interface  does  not  dissipate.  A  very  thin 
film  of  liquid  is  produced  between  the  projectile  and  the  target  plate. 

The  coefficient  of  friction  between  the  two  bodies  is  therefore  very  small 
so  the  frictional  forces  can  be  neglected.  This  conclusion  has  been  made 
by  Krafft  [12]  and  by  other  investigators. 

The  equation  of  motion  for  the  perforation  process  in  the  direction  of 
motion  is 

ft  (mV)  *  -P  Cl) 

where  F  is  the  resultant  force  acting  on  the  projectile,  m  is  the  instanta¬ 
neous  mass  of  the  projectile,  and  V  is  the  instantaneous  velocity  of  the 
projectile  in  the  medium. 

In  general,  the  force  F  is  the  resultant  of  three  main  components;  F^ 
the  inertial  force  of  the  target  material,  Fc  the  compressive  force,  and 
Fs  the  shearing  force-  Therefore, 


(2) 
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a?  <"v>  *  -  (Fi*Fc'Fs) 

The  particle  velocity  due  to  elastic  stress  wave  propagation  in  vhe 
target  is  neglected  in  this  analysis.  This  particle  velocity,  for  example, 
would  be  of  the  order  of  1  m/sec  for  a  0.22  inch  caliber  lead  bullet  moving 
at  a  velocity  of  400  m/sec,  e.g.  Davies  (13).  The  effect  of  the  propagation 
of  the  plastic  stress  vave  is  also  neglected  since  it  can  be  assumed  that 
its  velocity  would  be  smaller  than  the  velocity  of  the  projectile  for  most 
target  materials. 

The  projectile  therefore  transmits  kinetic  energy  only  to  the  mass  it 
displaces.  It  is  assumed  that  the  element  of  mass  of  the  target  material 
in  contact  with  the  projectile  is  set  into  motion  while  the  remainder  of 
the  target  material  remains  at  rest.  Each  mass  element  is  considered  to 
move  normal  to  the  surface  of  the  nose  of  the  projectile.  This  is  possible 
on  the  basis  that  the  target  material  is  compressible. 

Another  related  assumption  is  that  the  effective  mass  of  the  projectile 
increases  during  the  penetration  process  due  to  addition  of  the  target 
material  displaced  in  the  direction  of  motion.  Part  of  the  kinetic  energy 
imparted  to  the  added  mass  by  the  projectile  remains  stored  in  the  combined 
effective  mass  of  the  projectile  while  the  remainder  is  converted  to  plastic 
deformation  and  heat. 

In  the  first  stage  cf  the  penetration  process,  only  the  inertial 
and  compressive  forces  are  considered  to  act  on  the  projectile.  The  iner¬ 
tial  force  is  not  distributed  uniformly  on  the  projectile's  nose  surface 
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but  depends  on  the  shape  of  the  nose.  Equating  the  work  done  by  the 
reaction  of  the  inertial  force  on  the  target  material  to  the  change  of 
the  kinetic  energy  of  the  displaced  material  (neglecting  the  work  done  in 
changing  the  volume  of  a  mass  element  due  to  its  compressibility)  leads  to 


dFi»d\,  •  J  J"n2 


(3) 


where  dF,  is  the  normal  inertial  force  acting  on  an  element  area  dA  of 
in  n 

the  projectile,  dxn  is  the  displacement  of  a  mass  element  of  the  target 

material  normal  to  the  projectile  surface,  and  is  the  velocity  of  a 

mass  element  in  a  direction  normal  to  the  projectile  surface. 

According  to  the  above  assumptions,  the  mass  element  dm  of  the  target 

material  of  density  p  which  is  displaced  by  the  projectile  as  it  advances 

dx  would  be 
n 


dm  »  pdx  dA  (4) 

n  n 

Substituting  Eq.  (4)  into  Eq.  (3)  leads  to 

dFln  *  7  D<dVV„2  ‘5> 

can  be  determined  for  each  projectile's  shape  by  integrating  the 
component  of  dF^n  in  the  direction  of  motion  over  the  surface  of  the  pro¬ 
jectile. 


Nose  Shape 


Sphere 

Cone  (semi  apex  angle  a) 


.  2 

sin  a 


The  compressive  force  is  distributed  uniformly  in  the  direction  of 
>tion  on  the  projectile's  r.ose  surface  and  can  be  expressed  as  ^C“°CA- 
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A  is  the  projected  area  of  the  projectile's  nose  on  the  target  plate  and 
oc  is  the  ultimate  compressive  strength  of  the  target  plate  at  the  applied 
straining  rate.  The  equation  of  motion  for  the  first  stage,  Eq.  (2) ,  could 
therefore  be  written  as  follows: 


k  ("V) 


v  dm 
V  dt 


dV 

dt 


~  KpA.V2  -  o  A 
2  1  Cl 


(8) 


where  the  projected  area  A^  includes  the  effect  of  possible  flattening  of 
the  projectile's  nose.  A^  is  taken  to  be  constant  during  stage  1  and  is 
measured  by  the  entrance  opening.  In  the  actual  physical  process,  A^  would 
be  a  function  of  x  since  the  projectile  deforms  with  penetration  distance. 
This  effect  can  be  considered  in  an  overall  average  manner  to  within  the 
level  of  accuracy  of  this  analysis. 

The  inertial  force  F^  and  compressive  force  F£  are  considered  to  act 
on  the  effective  mass  of  the  projectile  which  includes  the  target  material 
displaced  by  the  projectile  and  moving  with  it  (Fig.  la).  The  problem  of 
physically  locating  the  added  mass  with  respect  to  the  projectile  is  not 
readily  resolved  since  compressibility  effects  would  have  a  large  influence. 
However,  this  determination  is  not  required  for  the  purpose  of  the  present 
analysis.  It  is  important  to  note  that  the  measure  of  the  penetration 
depth  of  the  projectile,  x,  is  that  of  the  combined  mass  and  not  that  of 
the  original  projectile  by  itself.  That  is,  the  value  of  x  is  the  distance 
from  the  initial  impact  surface  to  the  front  of  the  target  material  that 
is  moving  at  the  projectile  velocity  (Fig.  la).  The  effective  mass  of  the 
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combined  projectile  for  use  in  the  equation  of  motion  is  therefore 
m«*pAjX  whore  m,  is  the  original  mass  of  the  projectile. 

The  rate  of  change  of  the  effective  mass  of  the  projectile  would  be 

3T  ■  3F  *  <V  <9> 

Substituting  Eq.  (9)  and  the  relation 


dV  dV  dx  dV 

37  *  37  37  =  V  37 


into  Eq.  (8)  leads  to: 


pAjV2  *  (m0*pArx)  V  ~  ■  -  j  KpAjV2  -  ocA1  (10) 

Equation  (10)  can  be  solved  by  separation  of  variables  to  give 

„  .  .  2*K  1/2 

.2  °c  ^mo/pA1  .  o 

Vj(x)  *  '  (Vj  *p(1#0.5K)^  lViin,/pAl*x)  "  f(l*0.5K)  .  fll) 


where  V^(x)  is  the  velocity  of  the  combined  projtetile  and  added  mass. 

The  time  for  the  combined  projectile  to  penetrate  a  distance  x  can  bo 
calculated  by  numerical  integration  of  the  expression 


*(T  r dx 


(i 


/DA.  \ 


2.K 


fv  i  _ c _  ,  |  \  \  c 

1  i  ’p(l.O'.iK)  J  ym./pAj**;  *  P ( 1*0 . SKY 


-1/2 

dx 


(12) 
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The  forco-time  curve  for  the  first  stage  of  penetration  can  be  obtained 
from  Eqs.  (8),  (11)  and  (12).  The  first  stage  ends  when  x«h>b  (t*t^) 
and  the  process  of  shearing  starts.  The  value  of  the  plug  thickness  b,  which 
is  an  essential  factor  for  the  determination  of  the  various  quantities  at  the 
end  of  the  first  stage,  can  only  be  obtained  empiiically  at  the  present  stare 
of  development  of  the  analysis.  It  is  shown  in  (15)  that  the  ratio  b/h  is 
essentially  constant  for  a  given  projectile  and  target  material  within  the 
range  of  ordnance  velocities. 

In  the  second  stage  of  incipient  plugging,  the  inertial  force  continues 
to  act  on  the  projectile  and  can  be  expressed  as 

Fi2  "  J  KdA2V^ 

where  A.,  is  the  cross  sectional  area  of  the  cavity  in  the  second  stage. 

This  area  can,  in  general,  be  considered  a  function  of  x,  m  -^(x),  *n 
the  force  expressions  and  in  the  equations  of  motion.  The  experimental 
results  indicate  that  for  most  cases  A,  is  clcse  to  A^  so  that  the  complete 
cavity  can  be  considered  to  be  cylindrical  In  some  cases  there  is  apprecia¬ 
ble  enlargement,  i.e.  the  exit  diameter  is  much  larger  than  the  entrance 
diameter.  The  diameter  can  then  be  considered  to  vary  linearly  with  x  from 
l/j  at  .x«h  -b  to  at  x»h  to  give  a  quadratic  function  for  A^Cx).  The 
following  development  of  the  equations  of  motion  will  be  restricted  to  constant 
A,.  In  the  subsequent  comparisons  with  experiments,  [15],  A ^  and  were 
taken  to  be  equal  and  the  aveiage  of  the  entering  and  exit  areas.  Check 
studies  taking  A  =A  lx)  and  A  -A0(x)  showed  little  difference  m  the  final 

1  i  /  ** 
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results. 

The  shape  constant  K  in  the  inertial  force  expression  was  set  equal 
to  0.5  for  the  second  stage  since  for  standard  ordnance  projectiles  (not 
armour  piercing)  the  projectile's  nose  deforms  and  tends  towaid  a  spherical 
shape.  This  result  can  be  seen  from  ballistic  photographs  as  well  as  from 
the  geometry  of  the  ejected  plugs  and  the  deformed  projectiles.  This  change 
of  K  would  imply  an  artificial  discontinuity  in  the  force.  In  practice,  the 
discontinuity  is  very  small  and  not  observed  unless  the  time  increments  in 
the  computational  procedure  are  taken  to  be  very  small. 

The  compressive  force  alsc  acts  during  the  second  stage  of  penetration 
and  its  initial  value  would  be  F  «0  A_.  The  second  stage  ends  when  the  mass 
element  the  rear  side  of  the  target  plate  moves  at  the  same  velocity  of 
the  combined  projectile  and  effective  added  mass,  i.e.  x«h,  the  plate  thick¬ 
ness.  At  that  time  the  entire  target  material  forward  of  the  projectile 
moves  together  with  it  at  the  same  velocity.  Tne  force  F£  therefore  becomes 
zero  at  the  end  of  the  second  stage.  A  parabolic  function  for  F£2(x)  that 
meets  the  limiting  conditions  has  been  used  in  the  analysis, 

2 

Fc2(x)  »  ocA2  { 1  - ]  }  h-b  <  x  <  h  (13) 

The  alternative  choices  of  linear  functions  or  similar  forms  that  repre¬ 
sent  the  limiting  conditions  were  found  to  have  small  effect  on  the  calcu¬ 
lated  residual  velocities.  It  is  interesting  to  note  that  Eq.  (13)  does 
lead  to  force-time  curves  which  are  very  similar  to  those  obtained  in  the 
case  of  dynamic  punching,  e.g.  Dowling,  Harding  and  Campbell  [14]. 
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The  second  stage  of  penetration  is  also  characterized  by  a  shearing 
force.  This  acts  along  the  surface  of  that  part  of  the  plug  which  15 
moving  together  with  the  projectile,  i.e.  along  the  surface  ttD^ {x-(h-b) J 
where  D2  is  the  diameter  of  the  cavity  in  the  second  stage.  The  shear 
force  is  then  given  by 

Fs2(x)  *  tftDj  (x-(h-b)]  h-b  <  x  <  h  (14} 

The  shear  str.-.gth  of  metals  T  has  been  found  to  have  a  viscous  dependence 
on  strain  rate  at  very  high  rates  of  straining,  e.g.  [)<>],  (17] .  The  shear 
strength  can  be  taken  to  be  in  the  Bingham  form 

t  «  t.  ♦  UY  (15) 

where  p  is  the  coefficient  of  viscosity  and  y  is  the  shear  strain  rate. 

The  latter  cun  be  taken  as  V/e  where  e  is  the  radial  transition  distance 
between  the  plug  »nd  the  undeformed  target  material,  i.e.  the  width  of 
the  shear  zone  The  quantity  e  is  referred  to  as  the  "radial  clearance" 
in  dynamic  punching  problems  and  it  is  essentially  a  property  only  of  the 
target  material  at  high  rates  of  deformation  and  can  be  readily  obtained 
experimentally.  Analytical  expressions  for  e  can  be  deduced  from  the 
results  given  in  ( 1 8 |  and  (19|. 

The  equation  of  motion  for  the  second  stage  is  the  same  as  that 
derived  for  the  earlier  one  with  the  change  in  the  expression  for  the 
compressive  force  (13)  and  the  addition  of  the  shearing  force  (14),  (15). 
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The  equation  then  becomes,  for  the  case  of  constant  D2  and 


d 

a? 


(mV) 


1 

7 


KpAjV 


2 

(T.*U^1TD2[x-(h-b)J  -  <tcA2U-[^^-]  }  (16) 


where  h-b  <  x  <  h.  Substituting  Eq.  (9)  into  Eq.  (16)  leads  to 


dV2(x) 
Tk - 


pirD 

SK)DA2V  -  T.7TP2x  -  — 


UirD  (h-b) 

— Vx  ♦  - 1 - v 


[-  (1*0. 

♦  T«TTD2(h-b)  -  ocA2  (1  )  }  y/[  mj+pA2 


x)V] 


(17) 


where  m^  a  m«+pA, (h-b)  is  the  effective  mass  at  the  end  of  the  first  stage. 

The  time  t  for  the  combined  projectile  and  added  mass  to  reach  the  rear 
surface  of  the  turget  plate,  i.e.  x*h,  is  calculated  by  numerical  integration 
of 

x«h 

t2  »  j  ?  dx  (18) 

x*h-b 

V  can  be  calculated  numerically  by  a  computer  subroutine  from  Eq.  (17) . 

The  force  time  curve  can  be  obtained  from  Eqs.  (16),  (17)  and  (18). 

The  third  stage  cogences  when  the  entire  section  of  target  material 
forward  of  the  projectile  moves  together  with  it  as  a  rigid  body.  The 
effective  mass  is  then  m2*»m««pAh  where  A  is  the  average  cross  sectional  area 
of  the  entire  cavity.  During  stage  3  the  only  active  force  is  that  due  to 
the  shear  stresses  acting  over  the  surface  of  the  plug.  These  shearing 
stresses  are  considered  to  act  in  a  shear  lone  of  depth  e  around  the  plug. 
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The  displacement,  £,  of  the  combined  projectile  and  added  mass  system  with 
respect  to  the  plate  is  therefore  related  to  e  by 

5  -  Y®  09) 


where  y  is  the  shear  strain  in  the  effected  zone,  The  displacement  for 
material  failure,  E..,  is  reached  at  the  maximum  shear  strain  of  the  material, 
y^,  i.e.  £f.Yfe,  beyond  which  no  further  resisting  forces  act  on  the  moving 
system.  The  shear  strain  that  was  developed  in  the  second  stage  is  small  and 
could  be  neglected.  The  equation  of  motion  for  the  third  stage  is  therefore 


d2J. 

dt2 


F.  - 

a 


t  <■*  »<%  \ 


where  A^-Trl^b  and  is  the  average  cavity  diameter  in  the  second  stage. 
Using  Lq.  (15)  for  i,  Eq.  (20)  becomes 


r 
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The  force  during 

(24) 

The  time  duration  of  the  third  stage  t^  is  determined  by  the  time 
•required  for  the  displacement  £  to  reach  £^.  The  corresponding  velocity  at 
this  time  Vj  is  the  final  velocity  of  the  projectile  Y'^.  The  force-time 
relations  for  this  stage  can  then  be  determined  from  the  preceding  equations. 

The  total  time  for  the  perforation  process  is  the  sum  of  those  of  the 
three  stages  plus  the  time  required  for  the  plug  to  leave  the  target  plate. 
That  is 


wheie  ^  is  the  velocity  at  the  end  of  the  second  stage, 
this  stage  can  therefore  be  expressed  as 


,(t.  .  v  r)  [«P  (- 


t 


f 


*2  * 


(25) 


It  is  noted  that  this  time  would  correspond  to  that  of  full  ejection  of 
the  plug.  This  would  then  be  followed  by  ejection  of  fragments  correspond 
ing  to  the  effective  mass  added  during  stage  1  and  then  by  the  projectile 


itself. 
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Discussion 


r 


The  preceding  expressions  enable  the  post  perforation  velocity,  force- 
time  history,  and  contact  time  to  be  calculated  for  penetration  processes 
that  include  dishing,  plug  formation,  and  ductile  cavity  enlargement.  The 
relative  importance  of  the  mechanisms  considered  in  the  analysis  would  be 
determined  by  the  various  physical,  mechanical,  and  geometrical  parameters 
appearing  in  the  equations.  Of  these,  a  few  have  to  be  determined  empirically, 
namely  the  entrance  and  exit  hole  diameters  and  Dj  and  the  plug  length 
b.  Values  for  the  coefficient  of  viscosity  y  and  the  width  of  the  shear  tone 
e  can  be  obtained  from  the  results  of  other  investigations  and  modified  to 
suit  the  particular  ballistic  test  conditions. 

The  geometrical  measurements  to  be  taken  on  experimental  target  plates 

are  and  b.  Fairly  good  results  can  be  obtained  by  simply  setting 

°2  *  Davg  *  (D^Dj)  and  setting  for  the  area  corresponding  to 

Davg’  ^or  cascs  where  >>  ,  it  would  be  more  exact  to  calculate  A^  on 

the  basis  of  Dj ,  and  to  take  D^,  which  is  the  average  diameter  of  the  plug, 

to  be  the  linear  average  of  ind  over  the  plug  distance  b.  A  large 

number  of  measurements  on  a  variety  of  perforated  plates  [10],  [15]  has 

shown  that  the  ratios  D  /h,  and  b/h  are  essentially  constant  for  a  given 

avg  J 

target  piate  material  and  projectile  over  the  range  of  velocities  of  inter¬ 
est.  This  means  that  those  parameters  could  be  obtained  from  relatively 
lew  tests  and  then  used  in  the  equations  to  obtain  results  for  other  test 


conditions. 
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The  viscosity  of  materials  at  high  shearing  rates  has  been  determined 
by  various  rapid  loading  experiments,  e.g.  [14],  [16],  and  [17].  These 
results  could  be  used  for  u  in  the  analysis.  The  perforation  process  itself 
is  a  rapid  loading  experiment  so  u  could  be  considered  an  experimentally 
determined  material  property  within  the  framework  of  the  analysis.  That  is, 

V  for  a  given  plate  material  could  be  set  so  that  the  computed  results  would 
best  fit  the  test  results.  In  practice,  the  values  of  u  deduced  from  the 
ballistic  tests  and  those  that  had  been  obtained  from  more  direct  measurements 
are  generally  in  good  agreement  [IS). 

The  width  of  the  shear  zone,  e,  could  be  obtained  experimentally  by 
examination  of  etched  specimens  or  can  be  deduced  from  the  analyses  of  [18] 
and  [19i.  Again  both  methods  seemed  to  be  in  reasonable  agreement.  Either 
could  be  used  since  the  results  are  not  sensitive  to  the  exact  value  of  e. 

The  dependence  of  the  residual  velocity  on  the  properties  of  the  pro¬ 
jectile  and  target  plate  and  on  the  test  conditions  appear  to  be  in  confor¬ 
mity  with  general  observations  on  projectile  perforation.  A  typical  set  of 
f'  rcc-tirae,  velocity-time,  and  displacement-time  diagrams  obtained  from  the 
analysis  is  shown  in  Fig.  2.  This  example  was  calculated  for  the  case  of  a 
0.22  inch  caliber  lead  bullet  perforating  a  5.0  mm  thick  aluminum  alloy  plate. 
The  component  forces  throughout  the  three  stages  are  shown  in  Fig.  3  for  the 
same  case. 

The  absence  of  an  initial  rise  time  for  the  force  is  due  to  the  neglect 
of  the  shape  of  the  nose  of  the  projectile  on  the  rise  time  of  action  of  the 
compressive  and  inertial  forces.  Those  forces  are  assumed  to  act  immediately 
on  the  full  cross  section.  Consideration  of  this  effect  would  have  only  a 
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very  small  influence  on  the  overall  results.  The  decay  of  the  total  force 
during  the  first  stage  is  due  to  the  decrease  of  the  inertial  force.  The 
shear  force  acting  in  the  second  and  third  stages  is  seen  to  be  important 
and  dominant  in  the  last  part  of  the  perforation  process.  A  discontinuity 
in  the  force  (dotted  lines  in  Figs.  253)  would  appear  at  the  onset  of  the  third 
stage  due  to  the  removal  of  the  inertial  force.  Force  continuity  could  be 
maintained  ty  suitably  changing  the  value  of  e,  the  width  of  the  shear  rone. 
Detailed  examination  of  ejected  plugs  has  indicated  that  corresponding 
changes  in  e  do,  in  fact,  take  place  during  the  ejection  stage.  The  alter¬ 
ation  of  e  in  the  analysis  to  ensure  force  continuity  therefore  seems  to 
have  a  physical  basis. 

A  comparison  of  the  predicted  force-time  relations  with  those  obtained 
in  dynamic  punching  experiments,  e.g.  Fig.  4  from  [14],  shows  reasonable 
agreement.  In  those  curves,  the  time  over  which  the  force  decreases  after 
the  first  peak  would  correspond  to  stage  1,  i.e.  before  the  onset  of  plugging. 

There  are  important  differences  between  the  force-time  histories  obtain¬ 
ed  from  the  present  analysis  and  those  obtained  in  [11]  on  the  basis  of  an 
assumed  deceleration-time  function  of  the  projectile.  In  [11]  the  force 
tends  to  zero  in  the  last  part  of  the  process  while  the  present  results 
indicate  that  the  final  force  is  still  close  to  its  maximum  value. 

There  are  experimental  indications  that  the  force  is  large  at  the  last 
stage  of  perforation.  This  is  the  observation  made  in  [11]  and  other  bal¬ 
listic  experiments  that  the  velocity  drop  i.e.  the  differencj  between  initial 
and  final  velocities,  diminishes  for  velocities  slightly  in  excess  of  the 
ballistic  limit,  e.g.  Fig.  21  of  [11].  This  result  corresponds  to  the  obser- 
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vation  that  a  projectil'j  that  perforates  a  plate  under  conditions  slightly 
exceeding  the  ballistic  limit,  i.e.  a  higher  initial  velocity  or  thinner 
target  plate,  would  have  a  relatively  high  terminal  velocity,  e.g.  about 
20%  of  the  initial  value.  The  impulse  associated  with  the  end  process  of 
perforation  is  therefore  significant.  This  in  turn  implies  that  the  force 
acting  during  the  end  process  is  high  which  is  in  accordance  with  the  results 
shown  in  Figs.  2  and  3.  Another  consideration  is  that  the  time  interval  for 
the  third  stage,  tj,  is  found  to  be  sensitive  near  the  ballistic  limit,  and 
small  increases  in  thickness  or  decreases  of  velocity  would  inc* *a$e  t^  by 
a  relatively  large  amount.  This  effect  and  a  large  terminal  force  are  the 
apparent  causes  of  the  observations  on  the  velocity  drop  effect.  The  pre¬ 
dictions  of  terminal  velocities  based  on  the  present  analysis  (15]  do,  in 
fact,  show  that  the  velocity  drop  initially  decreases  with  increasing  initial 
velocity  as  indicated  in  Fig.  21  of  (11]. 
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Conclusions 

The  analysis  of  the  ballistic  perforation  problem  that  has  been 
developed  seems  capable  of  predicting  post  perforation  velocities,  contact 
times,  and  force-time  histories.  The  analysis  relies  on  certain  geometrical 
parameters  which  must  be  determined  empirically  at  this  stage  of  development. 
Certain  material  properties  are  not  well  established  and  could  be  determined 
by  extrapolation  of  other  results  or  by  experimental  observations.  Both 
these  and  the  empirical  geometrical  parameters  could  be  obtained  from  a  small 
number  of  tests.  Once  these  are  determined  for  a  given  projectile  and  target 
material,  predictions  can  be  obtained  over  a  wide  range  of  projectile  velo¬ 
cities  and  target  thicknesses.  Further  development  of  the  analysis  would 
be  to  determine  the  empirical  factors  by  basic  considerations. 
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Pig'  1  -  Schematics  of  the  stages  of  the  perforation  process 


Fig.  2  -  Calculated  example  of  displacement,  velocity,  and 

force  histories  for  the  three  stages  of  perforation 


Fig.  3  -  Force-time  relation  for  each  of  the  force 

components  for  the  three  stages  of  perforation 
(same  conditions  as  Fig.  2). 


Load/Displacement  curves  for  aluminum 
(from  Ref.  [14]).  '  r 


Fig.  4 
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